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ABSTRACT 



The background of microstructure and its effect on sound propagation 
is reviewed. The expendable bathythermograph (XBT) system and current 
digitizing procedures are evaluated to determine their adequacy in 
investigating vertical thermal microstructure. Several methods of 
microthermal analysis are proposed. Temperature gradient versus depth 
plots are used to analyze data from several ocean regions - Atlantic 
Coastal, Andaman Sea, Eastern North Pacific, and the northern boundary 
of the Gulf Stream. Analysis of these data suggest a water mass depend- 
ence of the thermal microstructure. The significance of the small-scale 
structure to Naval operations is discussed. 
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I . SIGNIFICANCE OF MICROTHERMAL STRUCTURE 



Naval interest in the thermal structure of the ocean is presently 
centered on gross features such as sea surface temperature, mixed layer 
depth, thermocline gradient, and depth to the bottom of the thermocline. 
Some attention has been given, especially in the Canadian Navy, to the 
transients in the mixed layer [Tully, 1964] and their possible associ- 
ation with internal waves. The above gross features are used primarily 
as a means of classifying thermal structure. 

Technological developments have led to measurements which indicate 
the thermal structure is made-up of layers of varying thickness; from 
less than a centimeter to tens of meters [Stommel and Federov, 1967; 
Osborn, 1969; for example]. 

Of the variables affecting sound velocity in the ocean, temperature 
has the greatest effect. The presence of variations in temperature, or 
temperature gradients, causes sound velocity fluctuations to occur. 
Research on the effects of the thermal layers on sound in the ocean show 
that further investigation in theory and experiment is required. Small 
scale sound velocity variations cause irregular bending of sound rays. 
The microstructure can affect attenuation, forward and backward scatter- 
ing, reverberation, reflection, and refraction which can converge or 
diverge the energy in a sound beam [Meyer and Romberg, 1963]. The 
physical features of the microthermal structure must be adequately 
described before their effects upon acoustic propagation can be fully 
understood . 
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The purpose of this thesis is to review the present knowledge con- 
cerning microstructure as measured by expendable bathythermographs (XBT) , 
to evaluate the capability of the XBT in measuring microthermal structure, 
to develop a scheme of quantitatively examining the microthermal structure, 
and to test this scheme on real data. 
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II. REVIEW OF PRESENT KNOWLEDGE 



A. MICROTHERMAL STRUCTURE IN THE OCEAN 

Prior to 1937, vertical thermal structure was measured at discrete 
points using reversing thermometers. This method provided a point-to- 
point analysis of the vertical structure at one location. Even using 
point-to-point measurements, it was clear that the temperature vari- 
ability decreased with depth in most regions. In fact, the most common 
analysis of the temperature and salinity data, T-S correlation, usually 
ignored the upper 100 meters or so. However, this same analysis showed 
that in many ocean regions, in deep water, (e.g., below 1000 to 2000 m) , 
the variability was on the order of the instrument error [Sverdrup, 
Johnson and Flemming, 1942]. This result was not unexpected because 
the gradients in the deep ocean became smaller. The discrete measure- 
ments of the Nansen cast are not suitable for studying small scale 
variability . 

The development of the mechanical bathythermograph (MBT) by Spilhaus 
[1938] provided the first convenient means of measuring a continuous 
profile of temperature. The mechanical bathythermograph was widely used 
during World War II. A large number of ba thythermograms with double 
traces, which could not be explained as instrument errors, were dis- 
covered [Urick and Searfoss , 1948]. The double traces indicated that 
small-scale horizontal and/or vertical temperature variations occurred 
between the lowering and raising of the MBT. Urick and Searfoss [1949] 
recognized that this small-scale microstructure could be significant to 
sound propagation, depending on the size and shape of the inhomogeneity. 
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Horizontal microthermal structure was measured by Holder in 1944 
using a thermopile mounted on a submarine. He detected horizontal 
temperature differences of 0.02F over distances of about 10 m [Smith, 
1967]. Measurements of these horizontal temperature variations were 
also made by Urick and Searfoss [1948, 1949] in July of 1948 near Key 
Wfest, Florida. With a thermocouple mounted on the conning tower and, 
later, its bow, a submarine proceeded at a constant speed and depth, 
while the probe measured the water temperature to a resolution of 
0.003C. 




depth- 25 H 





40 YO-^ 



Figure 1. Microthermal Variations at Various Depths 
[Urick and Searfoss, 1948]. 

Urick and Searfoss [1949] calculated the peak-to-peak root-mean- 
squara (rms) variation for peaks in excess of 20 percent of the total 
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variation in a record calling this the root-mean-square maximum (rmsm) . 

A significant difference between the value of rmsm in the mixed layer 
and in the thermocline was observed. In the mixed layer a mean rmsm 
value of 0.013C was observed. In the thermocline, the mean value was 
0.450C. The rmsm increased with depth in the mixed layer and decreased 
with depth in the thermocline. Applying Taylor's hypothesis, the mean 
horizontal thermal "patch size" was found to vary from 5 to 30 m in the 
mixed layer and from 30 to 200 m in the thermocline region [Urick and 
Sear foss , 1949 ] . 

Liebermann [1951] mounted a sensitive platinum-resistance thermom- 
eter and a thermocouple on the periscope of a submarine. The platinum- 
resistance element could resolve spatial temperature variations on the 
order of 10 cm in horizontal distances as the submarine proceeded at a 
constant speed of a few meters per second at a fixed depth (Figure 2). 

To characterize the temperature fluctuations, he used the autocorrelation 
function. "Patch size" was defined as the distance of the microthermal 
structure, within which the temperature variations maintained coherence. 

The function is given by: 



00 




( 1 ) 



0 



His results could be fit equally well by an exponential. 



R (p) = e 



( 2 ) 



or the Gaussian, 



R (p) - 



2 



( 3 ) 



where : 
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T (x) = temperatute at some point x; 

T (x+Pq) = temperature at another point displaced a distance 
from X ; 

p = autocorrelation distance; 

<p> = a = ’’patch size”. 

Whereas, the exponential function fitted the data points and decayed 
rapidly as p increased (Figure 2), it did not have a zero slope at the 
origin, thus, the exponential function was not appropriate for very 
small correlation distance (e.g., a few cm or less) because of heat 
conductivity in sea water. He proposed the Gaussian function to 
correct this difficulty, reasoning that thermal conductivity would 
prevent any discontinuity in the temperature field [Liebermann, 1951]. 

Liebermann assumed that the thermal microstructure represented a 
cloud of spheres with a Rayleigh distribution having a mean size value 
of the correlation distance and therefore yielding 

R (p) = R (a) = e’^ (4) 

Using this definition, the mean size of the inhomogeneities, a, was 
found to equal 60 cm in the ocean (Figure 2). 

Piip [1961] investigated the sound velocity microstructure in the 
sound channel off Bermuda with a NBS velocimeter. It was noted that a 
number of thin steps remained over periods of hours, some oscillating 
30 m in depth in a three-hour period. At the same station and same 
depth, sound velocity changes of one m/sec in a few hours were common 
and "nearly the rule”. The sound velocity axis was observed to change 
150 m in depth and 3.5 m/sec over a three-month period. Pipp showed 
that the base of the main thermocline (about 1000 to 1300 m) exhibited 
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relatively large and rapid changes in sound velocity in an inhomogeneous 
region in a constant state of fluctuation. 
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Figure 2. Temperature Fluctuations and Associated Correlation Function 

(After Liebermann) 

Small-scale temperature and salinity structure was measured with a 
salinity, temperature, and depth recorder (STD) by Stommel and Federov 
[1967] near Timor and Mindinao during July of 1965. Their measurements 
revealed that the vertical microstructure actually consisted of very 
thin high-gradient layers. The thin layers were separated by thicker 
layers of nearly homogeneous (isothermal or isohaline) water, i.e,, 
step-like in appearance on their STD traces. By following a discrete 
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feature from one cast to the next, homogeneous layers were found to be 
continuous from two to 20 km horizontally and to range in thickness from 
two to 40 m (Figure 3). A boot- like feature was found at the bottom of 
the mixed layer (Figure 4). At a depth of 200 to 300 m, the mean thermal 
gradient was about 0.040C/m. In the thin steps, the gradient was less 
than 0.003C/m, and in the homogeneous steps, the thermal gradient was 
greater than 0.350C/m. Stommel and Federov concluded that the ocean 
consisted of nearly homogeneous layers separated by steps with extremely 
high thermal gradients. 

Tate and Howe [1968], in the summer of 1966 in the Northwestern 
Atlantic, used an STD to look at the thermocline stratification in the 
Mediterranean Water intrusion at a depth of 1280 to 1500 m. They found 
layering with steps varying on the order of 0.17 to 0.35C in tempera- 
ture, from 0.02 to 0.055 o/oo in salinity, and from thickness of 15 to 
30 m (Figure 5). A step’s average temperature change of 0.25C produced 
a change of -0.044, while the average salinity step of 0.044 o/oo 
produced a change of +0.039. The resulting change of (-0.005) 
gave a slightly stable layer system in a situation which was clearly 
one of opposing salinity and temperature density gradients. 

Cooper and Stommel [1968], using an STD, found the main thermocline 
in the Sargasso Sea to consist of regular salinity and temperature steps 
of homogeneous layers three to five meters thick, separated by trans- 
itional layers 10 to 15 meters thick, with temperature changes of 0.30 
to 0.50C and salinity changes of 0.04 to 0.10 o/oo (Figure 6). The 
horizontal extent of these steps was between 400 to 1000 meters. The 
main thermocline contained up to 100 steps on each cast. 



14 





a. b. 



Figure 3. a. Photographic copy of actual STD trace (not retouched) of 

station 21 off Mindinao. Depth scale is 500 m, each small 
division being 5 m. The salinity curve, labeled S, has 
the smallest division being 0.05 o/oo. The temperature, 
labeled T, has a smallest division of 0.05C. The lamina 
B is so labeled. Because the pens cannot both traverse 
the same path mechanically, the temperature curve is 
offset downward by 5 m and always reads 5 m too deep. 

The salinity scale is not offset in depth. Uniform 
corrections for salinity of -0,9 o/oo and temperature of 
0.07C have not been applied. 

Figure 3.b. A composite chart-graph for lamina B off Mindinao. Each 
square is 0.05C wide (8.5 to 9.0C) and 50 m high. The 
small arrow indicates a depth of 400 m in each square. 
Numbers in the lower left hand corners are the station 
numbers. The distances from stations 18 to 25 and 26 to 
35 are about 10 nm. For lamina B, the large geographical 
extent of this small lamina is self-evident. 

[Stommel and Federov, 1967] 
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Figure 4. 



Figure 5. 






From the soundings in the mixed layer off Timor, gross dif 
ferences in the properties in different locations are self 
evident. There is considerable variation in vertical 
structure but each sounding shows a nearly homogeneous 
region near the surface (station 13) but some stations 

(like 7) have shallow irregular layers overlying the mixed 
layer. 



[Stommel and Federov, 1967] 

T&S 



T&S 




A typical step in Mediterranean Water [Tait and Howe, 1968] 
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Woods [1968, 1968a] investigated the microstructure in the summer 



thermocline off Malta using a temperature -gradient meter for obtaining 
the thermal gradient between two thermistors set 50, 25 or 10 cm apart 
on a vertical staff. The resolution of the thermal gradient was within 
O.OIC per thermistor separation distance. An additional thermistor 
recorded temperature. Woods plotted profiles of the thermal gradient 
and temperature versus depth (Figures 7, 8 and 9). Using time-lapse 
photography. Woods obtained pictures of streaks left in the wakes of 
free-dropping dye pellets for shear measurements. A dye-packet array 
tied between moored, submerged floats injecting dye into selected levels 
of the thermocline, indicated that the thermocline was divided into 
layers a few meters thick. These layers were characterized by weak 
thermal gradients, less than O.OOlC/m, and by turbulence with a mean 
velocity of less than one mm/sec. The layers were separated by sheets 
only centimeters thick. (Sometimes, he detected sheets up to a meter 
thick but they appeared as aggragates of the 10 cm thick sheets.) The 
sheets were characterized by strong thermal gradients up to 0.05C/cm 
and little or no turbulence. Mature sheets contained thin laminar-flow 
zones in their centers. Neighboring layers had dissimilar properties 
indicating different salinities, velocities, and turbidities. This 
suggested different trajectories and a very slow exchange of properties 
across the intervening sheets. If turbulence was small, flux exchange 
was small; therefore, the exchange had to rely on molecular processes, 
which were very slow indeed. 

Individual layers were followed horizontally for distances up to 
72 km by soundings from high-speed launches. Layers were identified on 
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Depth (metres) 




Figure 8. Temperature and Temperature -Gradient Soundings for Maltese 
Waters, 14 September 1967 [Woods, 1968a]. 

Temperature 0918-0929 local time, (l.e., GMT + 1). 

Temperature gradient 1132-1145 local time. 
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successive soundings for many hours, thus confirming the results of 
Stommel and Federov. 

The dye experiments showed that in a thermocline sheet there was a 
central laminar-flow region only a few centimeters thick. Dye was spread 
in this region by internal waves to form a ’’broad carpet of color”. This 
colored sheet made it possible to photograph the internal waves as they 
crossed it. Wave-lengths between 5 and 10 m were dominant and appeared 
in long, coherent trains, each associated with a single sheet. They had 
heights up to 1.0 m, a phase speed of 2.5 cm/sec and periods of about 
5 min. These waves only affected a single sheet, being much shorter than 
the long internal waves which move the entire thermocline. 

A sheet was defined as stable when it appeared smooth everywhere; 
this indicated laminar flow. When a wave passed over a sheet, shear had 
its greatest effects at the crest and trough. Wave shear and drift 
shear interacted either cancelling or reinforcing each other causing a 
sheet to become unstable and patchy. 

Wavelets, about 75 cm in length, caused banding to appear parallel 
to the crests and troughs. These wavelets generally formed, grew, and 
broke in a period of less than two minutes. Upon breaking in classical 
fashion, a second smaller breaker was generally thrown forward. The 
wavelets had a maximum height of about 20 cm. After breaking, a turbu- 
lence region remained for about 5 minutes, and the scar in the sheet 
remained for several hours. These scars were a common feature in 
thermocline sheets. The thinnest sheet observed to become unstable had 
a thickness of 3 cm; the thickest was 30 cm. The great majority were 
between 8 and 15 cm. 
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Transient thermoclines [Tally, 1964] which occur sporadically between 
the surface and the top of the thermocline were considered as further 
examples of sheets. Woods divides the thermocline into a half-dozen or 
so layers of low shear and moderate thermal gradient separated by thin 
liminar flow sheets on intense shear and thermal gradient. He discussed 
the important aspect of heat transfer through such a layered ocean. He 
suggested that the molecular heat transfer through the high gradient 
stable layers was negligible, but that^ heat transfer does take place 
through the turbulent scars from one layer to another as an effect of the 
breaking of the internal waves. 

The use of the Ramsey SVTP (sound velocity, temperature, and pressure) 
probe provides a digital output of these parameters. This instrument has 
an accuracy of + O.OIC and + 0.25 percent of depth with a precision of 
+ 0.08 m/sec in sound velocity [Lovett, 1968]. In view of the plus-or- 
minus-one count for electronic counters, small-scale features, i.e., 
less than + O.OIC and + 1 m in depth, are ambiguous. Lovett [1968] 
discussed the thermal gradients determined with this instrument (Figure 
10). The 10 m depth increments he used eliminated some microstructure; 
however, using 5 m increments introduced random readout errors. 

Ingham [1968] used an STD probe to investigate the mixed layer in 
tropical waters. The vertical salinity and temperature structures 
forming the "isothermal” and "isoha line" layers were generally analogous 
in the near-surface layers. The difference in the depths of these two 
layers exceeded 7.5 m in 45 percent of the casts and 17.5 m in 31 percent 
as shown by Figure 11 [Ingham, 1968]. Microstructure effects on sound 
velocity in the near-surface layers would depend on salinity in some ocean 
regions . 
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Figure 10 . Temperature Gradient per 10 m, Mean Gradient per 100 
Average Deviation from Mean versus Depth from Ramsey 
Data [Lovett, 1968 ]. 
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Figure 11 . Histogram of Depth Differences AZ of the Isothermal and 
Isohaline Layer Depths [Ingham, 1968 ]. 
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DEPTH IN METERS 



Pingree [1969] using an STD, recorded analogue signals on tape. He 
made fine-scale plots of salinity and temperature. In particular, he 
investigated the stability of the deep water, i.e., about 500 to 2000 m 
(Figure 12). 

Osborn [1969] using a free-falling temperature gradient meter 
obtained excellent profiles of temperature gradient (Figure 13). Off 
San Diego, the horizontal extent of the layers was greater than 750 
meters in the seasonal thermocline, but only a few hundred meters at 
depths greater than 400 meters. He noted the layers were, in fact, 
made up of even finer layers with a thickness of only a few centimeters. 
His average sheet was about one meter thick. 

Recent work by Neshbya, Neal, and Denner [1969] using the expendable 
bathythermograph manufactured by General Motors Defense Laboratories, 
obtained records of temperature steps (Figures 14, 15 and 16). While 
on Ice Station T-3 in the Arctic (near 84-38N, 128-22W) , they lowered 
XBT probes by hand. The Ice Station results are the first XBT micro- 
structure records without ship motion. 

B. APPLICATIONS 

Urick and Searfoss [1949] calculated theoretically that a sound ray 
would be bent through an angle of only 3.7 sec by the average inhomo- 
geneity in the mixed layer but as much as 2.1 min of arc in the 
thermocline . 

Sheehy [1950] measured sonar signal fluctuation as functions of depth 
and range. He found the coefficient of variation, defined as the stand- 
ard deviation of the amplitudes expressed as a percentage of the mean 
amplitude, varied as the square root of range and did not appear to be a 
function of depth. 
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Figure 12. T-S Curve from 1500 to 1800 m with a Value Every Meter. 

Circles Represent Increments of 2 m in Depth [Pingree, 
1969], 
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Figure 14. Vertical Profile of Temperature Under Ice Station T-3 
(84-38N, 128-22W) 19 March 1969. [Neshbya, Neal and 
Denner, 1969]. 
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Liebermann [1951] described the influence of the microstructure on 
acoustic propagation He wrote the sharpness of temperature boundaries 
determines influence of the microstructure on acoustic propagation. 

When temperature changes occur within a distance of less than one 
acoustic wavelength, the temperature boundaries may reflect sound 
energy depending upon the size and shape of the "patch" The tempera- 
ture microstructure may reflect, refract, scatter, or focus acoustic 
waves r The microthermal structure was considered to be made up of 
randomly sized "spheres", spaced randomly throughout the medium, with 
a scattering cross section dependent upon the autocorrelation of the 
temperature inhomogeneities, R (p) . Comparing derived scattering cross 
section with actual scattering in the ocean (ie., volume reverberation), 
Liebermann concluded the observed reverberation was higher than one 
would expect from thermal inhomogeneities He attributed volume 
reverberation to biological scatters. 

Liebermann found that refraction effects of the inhomogeneities 
resulted in warping of thp wave fronts, and caused the wave front to 
be a complex surface exhibiting both convergence and divergence. The 
temporal motion of the inhomogeneities resulted in sound "scintillation" 
at the receiver. 

Mintzer [1953] found Sheehy*s coefficient of variation to be depend- 
ent upon the three-halves power of the range rather than the one-half 
power of range. 

Skudrzyk [1957] discussed the focusing and defocusing effect within 

2 

the focusing range r^ given by r^ = ka where a is the "patch" radius 
and k is the acoustic wavenumber At ranges greater than r^ , in the 
interference range, phase interference results (Figures 17 andl8). 
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Figure 17. Illustration of the Various Types of Acoustic Scattering 
Due to a Thermal Patch in the Sea [Whitemarsh, Skudrzyk 
and Urick, 1957]. 




Figure 18. The Phase Angle Between the Incident and the Scattered 
Pressure and the Physical Nature of the Scattering as a 
Function of the Elange [Skudrzyk, 1957]. 
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Using this concept, Skudrzyk (Figure 19) shows the focusing range r^ as 
a function of frequency and patch size. (For example, a 50 kHz frequency 
sonar has a focusing range of 75 meters, if the patch size is assumed to 
be 60 centimeters.) For finding patch size, Skrudzyk used the Urick- 
Searfoss [1948] data to show patch size was twice the depth (Figure 20). 
Further, Skudrzyk started theoretical investigation with spheroidal 
patches, which he found to have a considerably stronger focusing effect 
than the spherical patches. 

Whitemarsh, Skudrzyk, and Urick [1957] did additional work with the 
spheriodal theory. The correlation functions previously used were dis- 
cussed. It was concluded that the present correlation functions were 
not suitable and either had to be replaced or redefined. The authors 
argued if the temperature microstructure was produced by turbulent 
motion in the sea, then the Kolmogorov Law for statistical equilibrium 
must hold. This, they show, was in excellent agreement with their data. 
Lending support to continued studies such as are presented in this 
thesis, their final conclusion states that the 

"... greatest gap seems to lie in the microthermal 
statistics of the ocean. If these were known better 
than they are now, apparently the amplitude statistics 
of the fluctuations of sound transmitted through the 
sea could be predicted roughly, but with some assurance."^ 

Works by Peterson [1965], Barakus [1968] and DiNapoli [1969] have 
considered the normal mode approach to the interference range solution 
for the temperature inhomogeneities and have also investigated the 
effects of internal waves on sound in stratified mediums. 



Whitemarsh, D. C., Skudrzyk, E., and Urick, R. J. , "Forward Scattering 
of Sound in the Sea and Its Correlation with the Temperature Micro- 
structure, J. Acous. Soc ■ Am. , v. 29, p. 1141, October 1957. 
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Figure 19. Focusing Range (meters) as a Function of Frequency and 
Patch Size (meters) [Skudrzky, 1957]. 
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Figure 20. Patch Diameter Versus Depth [after Skudrzyk, 1957]. 
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C. INADEQUACIES OF THE PRESENT THEORY 



Since the thermal microstructure is not related simply to the gross 
thermal features measured on MBTs , there is presently no easy way to 
predict, on a day-to-day basis, the effects of microthermal structure 
on the propagation of sound in the sea. The theories, so far, have 
developed in two classes. 

The first theories [Liebermann, 1951; Mintzer, 1953, 1953a, 1954; 
and Skudrzyk, 1957] concerned spherical and, more recently, spheriodal 
shaped inhomogeneities throughout the ocean. Part A above has dis- 
cussed the present knowledge of microstructure showing that the long 
thin sheets and layers do exist as microthermal structure. These first 
acoustical theories do not apply to what is now described as the micro- 
structure of the ocean. 

The theoretical treatment by Whitemarsh, Skudrzyk and Urick [1957], 
considering speriodal patches and Kolomogorov ’ s turbulence theories, 
could not find any real theoretical explanation of what happens in the 
sea when sound passes through an inhomogeneity. The investigations 
using normal modes have not produced any practical operational results. 

Lee [1961] using ray path theory, demonstrated the effect of a 
typical sinusiodal internal wave on sound intensity in a three- layered 
model. For the medium with an internal wave, intensity variations up 
to 22 db were found in distances of less than one internal wavelength, 
as shown by Figures 21 and 22. In contrast, the same medium, without 
an internal wave, produced a 5 db variation, including a 2.5 db 
spreading loss, over the same distance. 

Another study investigation the effects of sound velocity micro- 
structure was conducted by FNWC [Wolff, Tatro and Megehee , 1967]. A 



t 



35 




•itTAMCi IN fccr 



Figure 21. Diagram of Sound Rays through a Medium Which has an 
Internal Wave on the Thermocline [Lee, 1961]. 




Figure 22. Sound lavel in the Medium with an Internal Wave on the 

Thermocline. The db reference level is that corresponding 
to a sound level of 60 db at 1 foct from the directional 
•ource (X equals 0, 2 equals 10 ft) along the horlrontal. • 
The field is contoured at intervals of 2.5 db [lee, 1961). 
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sophisticated ray-tracing program [Ayers, Wolff, Carstensen and Ayers, 1966] 
was used for the analysis. Three complete ray trace runs were made to 
the first convergence zone using 120, 54, and 20 digitized sound velo- 
city points from one SVTP trace containing sound velocity microstructure. 

The 120-point sound velocity profile preserved all the '^wiggles" or 
microstructure, the 54-point profile used values at every 50 meter of 
depth, and the 20-point profile used the standard oceanographic depth 
values for sound velocity (Figure 23). The range to the first conver- 
gence zone agreed within 150 yards for all three cases. Based on the 
ray trace runs only, it was concluded the sound velocity microstructure 
was "essentially self-cancelling" and did not affect the convergence 
zone parameters [Wolff, Tatro, and Meghee , 1967]. Propagation loss and 
reverberation effects were not reported. 
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III. MEASURING MICROTHERMAL STRUCTURE 



A. THE BATHYTHERMOGRAPH 

1 . Mechanical Bathythermograph 

At the present time temperature is the most readily observed 
variable in the upper few hundred meters of the ocean. A major con- 
tributor to this situation is a singularly successful oceanographic 
instrument known as the mechanical bathythermograph or MBT. Subsequent 
to its introduction by Rossby and Montgomery [1934] and perfection by 
Spilhaus [1938], the MBT has become one of the basic instruments of 
almost all oceanographic agencies because of its simplicity of 
operation, reliability and relatively low cost. The result has been 
that a large number of bathythermograms have been collected. 

Researchers during World War II, using the MBT, observed micro- 
thermal structure on their traces. No theory had been formulated which 
took into account the effects caused by microthermal structure on under- 
water sound transmission. It was recognized, however, that temperature 
microstructure must be considered in any explanation of observed under- 
water sound transmission [National Defense Research Committee, 1946], 

2 . The Expendable Bathythermograph 
a. Capabilities 

Presently the U. S. Navy is using the XBT system manufactured 
by the Sippican Corporation of Marion, Massachusetts. This system con- 
sists of an expendable shipboard cannister and probe, shipboard launcher, 
and a strip chart recorder. Data so obtained are submitted to FNWC and 
were used in this thesis. 
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The original Bureau of Ships* specifications [Denner, 1966] 



for an XBT system are given in Table I. 

TABLE I 



Temperature range 
Temperature accuracy 
Depth range 
Depth accuracy 

Ship’s speed when launching 



28F to 96F 
+ 0.4F 
0 to 1000 ft 

+ or 15 ft, whichever 
is greater 

+ 63^ of a temperature 
step after 3 ft depth 
change 



The expendable Sippican probe currently used has a 
thermistor sensor in its nose. The probe is connected electrically 
to the recorder by two conductor wires which pay out from both the 
probe and the canister. The cannister in the launcher is connected 
to the recorder by a cable. There are no splices in the wire. The 
probe’s weight is a critical factor in determining its fall-rate. 

The Sippican system was found to be within + 3 gms both with and 
without its internal copper wire. (The deviation was from 0.3^ to 
0.64^ of the total weight of the probe.) The observed average fall- 
velocity of the probe was 20.3 + 0.2 ft/sec, A probe, with its 
internal wire, had an average terminal velocity of 19.9 + 0.7 ft/sec. 
Terminal velocity was reached at 18.0 + 0,5 ft [Gouzie, Sanders and 
Littlehale, 1966]. 

The present XBT strip chart recorder is a Mark 2A model 
using pressure sensitive paper for recording the thermal profile. The 
paper moves vertically on the recorder at a rate of 3.25 in/min. In one 
minute, the XBT, falling at 20 ft/sec, will have reached a depth of 1200 
ft. Thus, 1200 ft will be recorded on 3.25 inches of chart paper, or 
0.0027 of an inch of chart paper is used for every foot of fall. 
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The chart paper records horizontally from 28F to 96F in 
7.0 inches, or about 0.103 of an inch of chart paper per l.OF. Table 
II shows the resolution capabilities of the paper using different 
recording intervals. 



TABLE II 

Chart Spacing (inches) 

0.01 

0.005 

0.001 



Resolution 

Temperature (F) Depth (ft) 

0.0971 3.690 

0.0476 1.045 

0.0097 .369 



The response time for the thermistor is 110 milliseconds 
nominal and 130 milliseconds maximum [Choate, 1970 and Demeo, 1968]. 

A time constant-temperature change plot shows that 98.2 percent of a 
temperature change has been recorded in 10.4 feet of fall, 83 percent 
recorded in 5.2 feet of fall, 68 percent in two feet of fall and 35 
percent in one foot of fall (Figure 24). 

The thermistor, because of its finite size, has thermal 
inertia which causes it to act as a low pass filter. As a result, 
fine structure is lost or at best smoothed out. The original trace 
could be reconstructed if the filtering characteristics (of both phase 
and amplitude) of the thermistor were known. These could be obtained 
experimentally . 

b. Data Handling 

Present procedures for handling XBT data are indicated by 
the log sheets which accompany the XBT’s in the shipping container. 
First, a BATHY or BAXBT message is required to be submitted to the 
appropriate addressee using a PRIORITY precedence message [National 
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Figure 24. Temperature Change as Functions of Time, Time Constants, and Depth for Various 
Temperature S te ps . 



Oceanographic Data Center, 1968, 1968a]. Upon receipt of the BATHY or 
BAXBT message the receiving agency uses this as synoptic data for sub- 
jective ocean forecasts in his locality. The BATHY or BAXBT is for- 
warded through the Naval Environmental Data Network to FNWC in Monterey. 
FNWC uses the message reports to update its ocean fields for objective 
forecasting in various oceanographic and meteorological products. 

Upon completion of a cruise the ship which took the XBT 
cast mails the original trace and log sheet to FNWC for digitizing. 

The original XBT trace is digitized for use by FNWC, and, when 
requested, the results of the digitizing are provided to the XBT 
submitting activity. Upon completion of digitizing, the original XBT 
traces are forwarded to the National Oceanographic Data Center (NODC) 
for storage in their archives. NODC also receives a copy of the digital 
data on tape. 

Present procedures at FNWC use DIGITBT and several other 
specialized XBT data handling programs. New XBT digitized data are 
stacked on the XBT Master File (presently at a rate of about 500 traces 
per week) . 

B. OTHER SENSORS 

Sensors for measuring microthermal structure, other than those 
discussed elsewhere in this thesis, are numerous. 

From a fixed platform, like the NEL Oceanographic Tower, thermistor 
strings and arrays are used. Bottom devices, such as thermistors, 
attached to the shore or a tower by a wire, have been employed to give 
continuous temperature records. 



43 



From submarines, thermopiles, thermocouples, thermistors, platinum 
resistance thermometers, and thin film thermometers have been employed. 
Almost any oceanographic electronic or mechanical thermometer could be 
used on a submarine . 

C, DIGITIZING PROCEDURES 

At the present time, all XBT traces are submitted to Fleet Numerical 
Weather Central (FNWC) in Monterey, California, the only digitizing 
center for these strip chart analogue traces [National Oceanographic 
Data Center, 1968, 1968a; and Bauer, 1969]. 

At the Point Pinos Annex of FNWC, two Calma Model 480 flat-bed 
digital recorders are used to digitize XBT traces in 0.01 inch 
increments. These digital recorders are operated in a mode which makes 
the recorded intervals dependent on the tracing speed. The Calma 
output recorded on digital magnetic tape containing binary mode 
characters at a density of 556 bits per inch (BPI) with even parity 
and variable length records. 

FNWC program DIGITBT, written by Lieutenant N. L. Perkins, USN, 
in FORTRAN IV for FNWC's CDC 6500 computer, is used currently for 
the reduction of data from XBT traces. FNWC tape input, tape output, 
and some other subroutines of DIGITBT are written in COMPASS language. 

In order to analyze microthermal structure, DIGITBT was extensively 
modified. The new program, referred to as MICROXBT, is included as 

Appendix A. The basic procedures for DIGITBT and MICROXBT programs 
are the same . 

DIGITBT and MICROXBT programs both convert the X and Y measure- 
ments of the digitizer into depth (X) versus temperature (Y) values. 
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Corrections are made for nonlinearities produced by the XBT recorder in 
both the depth and temperature directions. All shipboard XBT recorders 
are assumed to produce error-free traces. The 62 . OF temperature line 
drawn on the trace when the probe is loaded is used as a baseline for 
digitizing. No correction is made if this baseline temperature line 
is not the same as the 62, OF line printed on the trace paper. An 
eye- interpolated value of the XBT-drawn baseline temperature on the 
printed temperature grid is noted as a part of the identification 
record . 

The rigid procedures for digitizing require that exactly one 
identification (ID) record and two digitized trace records for each 
XBT be recorded on the tape. The ID record information contains the 
data (day, month, year); time (hour and minute); latitude and longi- 
tude (degree, minute, hemisphere); baseline temperature (to nearest 
O.IC or O.IF depending on the grid type); printed-grid type; and maxi- 
mum depth for the probe (coded as 15 for 1500 ft XBT or 25 for the 
2500 ft XBT). Each trace record contains pairs of increments (AX, AY) 
read from the baseline to the true surface temperature, a point "flag” 
and then pairs of values (AX, AY) along the trace to a cut-off depth. 

Each trace is digitized twice. The complete procedure is described 
by Bauer [ 1969] . 

To process digitized XBT data using the DIGITBT program, the two 
trace records are compared point-by-point in a deviation test. If 
the X and Y paired values (i.e., depth and temperature) differ by more than 
0.035 inches, DIGITBT rejects the whole XBT record and the XBT trace must 
be completely re-digitized. If the XBT record passes the deviation test. 
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the two traces are averaged, point-by-point. The resultant trace is 
used for the rest of the program. 

In the following the term "inflection point" will be used to pre- 
serve consistency with the author of DIGITBT. Inflection points are 
those points which meet the programs^ criteria as significant points 
which are the programs output. 

Both DIGITBT and MICROXBT select the first and last observations 
as inflection points in the output list. The two programs differ at 
the inflection point choice stage. Other observations in DIGITBT are 
screened in groups, and, for each group, the single observation is 
found which is most different from a linear interpolation between the 
two end points of this group of observations. Each observation in a 
group is measured to find how out-of-line it is in temperature and in 
depth, and the smaller of these two is considered the out-of-line 
measure. All observations in the group are then compared to find the 
one with the largest out-of-line measure. If that value exceeds 
0.0150 for an inch, that observation is selected for output in DIGITBT; 
otherwise, the group expands to include the next value and the inflec- 
tion test is repeated until an inflection point is found. 

The selection of inflection points was modified in MICROXBT to use 
a minimum of 0.0001 of an inch for the out-of-line measure. This pro- 
duces an output containing every inflection point, which is nonlinear 
between two consecutive points in the digitized record. Any observations 
which cannot be reconstructed by linear interpolation are not required 
and are removed from the output listing in both DIGITBT and MICROXBT. 

DIGITBT converts depth units from digitized inches to ft of depth 
by linear interpolation using a depth table and from ft to m by formula. 
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Temperature units are converted from digitized inches to F by a tempera- 
ture table and from F to C by formula. These depth and temperature 
tables were constructed to account for the nonlinearities between the 
recorder and the actual trace. Values of depth are rounded to the 
nearest m and temperature to the nearest O.OIC in DIGITBT. MICROXBT 
has eliminated the ft to m and F to C conversions. Depths are rounded 
to the nearest 0.1 ft and temperature to the nearest O.OOIF. 

After the conversion from ft to the nearest m, DIGITBT removes, 
for a second time, any of the rounded-off values which are nonlinear 
between two neighboring points. If depth values appear as duplicates, 
one value is omitted in the output listing. In MICROXBT, there are no 
duplicate values because 0.1 ft in depth is well beyond the hundredth 
of an inch digitizing capability of the digitizer. Thus, by reading 
depths in tenths of a ft and temperature in thousandths of a degree F, 
MICROXBT looses no values from duplicate depths, and the possibility of 
introducing false linearities is reduced considerably. For these rea- 
sons, MICROXBT gives a much finer digitized record than DIGITBT. 

The DIGITBT output tape is written in binary for the GDC 6500; the 
MICROXBT output tape is written for the IBM 360/67 at the Naval Post- 
graduate School (NFS). Appendix C contains the job control cards 
(JCL) required for using a GDC 6500 output tape from MICROXBT with 
external BCD at a density of 556 BPI as in input to a separate micro- 
thermal analysis program at NFS. 



47 



IV. MODELS FOR EXAMINING MICRQTHERMAL STRUCTURE 



All models and methods used in this thesis required the output tape 
from MICROXBT to be used as a direct input on the seven channel tape 
unit. The XBT ID record information was read first and if the recorded 
XBT was from a desired geographical region the data (depth and tempera- 
ture) were analyzed for microthermal structure at NFS. 

A ”best*' value for the temperature structure was found by using the 
depth inflection points from MICROXBT. For this *’best*' array, tempera- 
ture was linearily interpolated to a value for every foot of depth. The 
mixed layer depth (MID) was determined, using inflection points from 
MICROXBT, as that depth immediately above the depth where the temperature 
was 2. OF different from that at the surface [Naval Weather Service, 1967]. 

Several microstructure models were tried in order to perform an 
objective analysis of the oceanic microthermal structure. The use of 
a numerical temperature derivative model was found to give the most 
satisfactory presentation of the microthermal structure. 

A. DERIVATIVE MODEL 

1 • Mathematical Model 

The most significant feature of the thermal structure is the 
thermal gradient defined as the derivative of the temperature with 
respect to depth (dT/dZ). By plotting the gradient versus depth (Z) , 
the thermal structure of the water column is emphasized. Furthermore, 
a parameter such as the thermal gradient is related to stability, 
entrophy generation, turbulent heat transport, and the Brunt-Vaisa la* 
frequency (Osborn, 1969]. 
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2 . Experimental Methods 



There are two methods for experimentally determining the thermal 
gradient presently in use. 

a. Electronic 

A resistance thermometer in an electronic bridge passing 
through water produces a signal which is proportional to temperature as 
a function of time. By differentiating this signal with respect to time 
(which is analogous to differentiating with respect to distance), the 
values of gradient versus distance the probe travels are found. 

This method was successfully employed by Grant, Moillet and 
Vogel [1968], in the horizontal direction for temperature and turbulence 
measurements. A thin-film resistance thermometer, mounted on the bow of 
a submarine, measured temperature versus horizontal distance traveled. 
The time derivative was taken to produce a direct plot of dT/dX, the 
horizontal thermal gradients. 

Osborn [1969] used a thermistor, capacitively coupled to 
two differentiating amplifiers, to measure the vertical temperature 
gradient with his free-falling gradiometers . His resulting output 
(Figure 13) shows a gradient versus depth record for the North Pacific 
off the San Diego Trough. 

b. Two-probe System 

Temperature gradients may be found by using two vertically 
spaced thermistors. The temperature difference is a measure of the mean 
gradient between them. Such a gradient-meter was used by Woods [1968, 
1968a] who obtained gradient at various depths (Figures 7, 8, and 9). 
This method is limited to the resolving power of the thermistors and is 
a function of thermistor spacing. 
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3. Numerical 



The numerical method is based on the derivative, dT/dZ, for the 
gradient but uses the finite difference approximation. The derivative 
of the temperature at any depth T(Z^) with respect to depth Z is defined 
as : 



dZ 



lim 

Z.-Z 



T(Zp - T(Z) 

zT - Z 
1 



(5) 



which can be shown to equal: 

^ ^ lim ^ 
dZ AZ-*0 AZ 

( 6 ) 

lim T(Z + AZ) - T(Z) 

AZ-0 AZ 

For this thesis, AZ = 1 foot, which yields the finite difference approxi- 
mation for the vertical temperature derivative: 



^ ^ AT. 
dZ 



(7) 



B. MISCELLANEOUS METHODS 

1- Statistical Bias and Variance 

For this method, a standard or smooth curve was required. 

Several standard curves were investigated. The variance and bias were 
calculated between the standard and the "best" curves. The standard was 
found by computing the temperature values for a series of specified depths. 

The first standard curve investigated used MLD and the FNWC 
significant levels (surface, 100, 200, 300, 400, 600, 800, and 1200 ft). 

The temperature values were found using the "best" temperature values at 
those depths. Next, an array of temperature values was found at foot 
increments by linear interpolation between the temperature values at the 
significant levels or the MLD. This was designated the FNWC significant 
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level curve. Calculations of bias and variance did not reveal any cor- 
relation with the microthermal structure easily visible on the XBT trace. 

The second standard curve investigated used the temperatures at 
the standard oceanographic station depths (surface, 10, 20, 30, 50, 70, 
100, 150, 200, 250, 300, 400, 500, and 800 m) . Since these depths are 
metric, a linear interpolation was done between the two neighboring 
values on the "best” curve. Values of temperature for every foot were 
found by linear interpolations between the standard depth temperature 
values. Again, no significant correlation was noted. 

The third standard curve used thse standard depth temperature 
values but used these standard depth values as an input to the LGINTP 
subroutine [Denner, 1969]. This subroutine uses the four standard depth 
points neighboring the depth in question to make an "average" curve from 
an upward and downward parabolic interpolation and a linear interpolation. 
This scheme produces a smoothing of points on the curve. This technique 
is usually used for interpolating standard depth oceanographic data 
[Rattray, 1965]. These values, after computing the bias and variance, 
again gave no correlation with the microthermal structure. 

The LGINTP subroutine was not used with the FNWC significant 
level values for smoothing since the FNWC fields assume linear inter- 
polation . 

Other standards considered, but not investigated, included a 
Gaussian standard, using the technique for approximating the temperature 
structure proposed by Grosfils [1968] and an exponential standard. 

2 . Gradient Distribution Method 

For investigating the temperature structure, the "best" tempera- 
ture data were separated into the FNWC significant levels. For each 
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significant level, the temperature change from one depth to the next 
were calculated and then sorted in numerical ascending order with the 
IBM subroutine SHSORT, With these values of temperature change per foot 
or gradient (F/ft), the temperature structure was now characterized by 
the individual thermal gradient values and a total temperature change AT. 
Because of the range of temperature values in the thermocline, this region 
of the XBT trace was chosen for developing this method. 

The capabilities of the XBT indicate a gradient interval of 0.02 
F/ft is easily obtained with a high degree of confidence [Arthur D. 

Little, 1966]. Using 0.02 F/ft for the class interval, the frequency 
distribution of gradients was plotted for various XBT traces by signifi- 
cant levels. As expected, the greatest range in gradient values usually 
occurred in the significant level just below the mixed layer where the 
range of temperature values AT was the greatest. 

This method appears to give a good indication of the variability 
of the thermal structure in the depth intervals investigated. Several 
examples of the different types of gradient distributions, using histo- 
grams and a cumulative frequency diagram, are shown in Figures 25 through 
31. The actual analyzed XBT is also shown on the left of each figure. 
Looking at data off San Diego, the first three gradient distribution 
diagrams (Figures 25, 26 and 27) show fairly similar cumulative fre- 
quency distributions. The next two (Figures 28 and 29) show nearly 
identical distributions for water inside the Gulf Stream. The next two 
gradient distributions (Figures 30 and 31) show water in the Gulf Stream's 
frontal region. These last two distributions are not identical for many 
reasons, but mainly because the positive gradient region in the second 
example is below the 100 to 200 foot interval shown in the gradient 
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Figure 31. Gradient Distribution for the 100 to 200 Foot Interval for Water in the Gulf Stream' 
North Wall. 



distribution plot. These results indicate that the small scale structure 
may be water mass related. However, this could not be fully examined in 
the time alloted for this thesis but is recommended for further 
investigation . 
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V . APPLICATION OF NUMERICAL DERIVATIVE MODEL 



This chapter describes the application of the numerical derivative 
model or method (NDM) to several cases of real and assumed data. 

A. BASIC CASES 

Graphs are used to show the features of the NDM. For each station, 
temperature versus depth (T-Z) and temperature gradient versus depth 
(G-Z) plots were drawn by a CalComp plotter using the computer program 
with the DRAW subroutine in Appendix B. The scales are constant for all 
plots. The grid squares were drawn for 100-ft depth increments from the 
surface to a depth of 900 ft. On the T-Z plots (reporduced digitized XBT 
data), temperature was scaled at lOF per grid line. On the G~Z plots, the 
gradient was scaled to O.lF/ft per grid interval, with the axis offset so 
that two-thirds of the plot was allocated to negative gradients. This 
scaling was chosen so that each hundredth of an inch, the resolving 
power of the CalComp plotter, represented one foot of depth and O.OIF 
of temperature or O.OOlF/ft of temperature gradient. The plots included 
here have been reduced in size, but the relative proportions remain the 
same . 

1 . Linear Gradients 

For a temperature distribution with linear thermal gradients, 
the G-Z plot is characterized by steps indicating a constant gradient 
during a temperature change. For this thesis the strength of a gradient 
is defined as the absolute value of the gradient (in F/f t units) and is 
described as strong when greater than 0.3F/ft or weak when less than 
0.05F/ft. The term signed strength is defined as the true gradient 
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value and will thus distinguish absolute and true gradient values. The 
magnitude of a gradient is defined as the vertical extent of the gradient 
measured in feet. The product of magnitude and signed strength gives the 
total temperature change for the extent of the gradient. 

a. Layer Structure 

Figure 32 shows a two- layer temperature structure model with 
a T-Z plot on the left and a G-Z plot on the right. The strength of the 
upper layer is 0.05F/ft with a magnitude of 200 ft. The lower layer is 
weaker (0.02F/ft compared to 0.05F/ft) than the upper step but the magni- 
tude is greater (600 ft to 200 ft). 

For the multi-layered model (Figure 33) several linear 
gradient steps are shown. In the gradient region between 70 and 90 ft, 
the gradient is strong and goes off scale (Figure 33b). 

b. Step Structure 

A step structure is characterized by thin strong gradient 
regions. The model (Figure 34) shows clearly that the four steps bet- 
ween 60 and 100 ft are stronger than any other steps shown. These four 
steps also have the smallest magnitude (five ft) in the figure. The 

steps at 200, 300, and 400 ft all have the same strength but increase 
in magnitude with depth. 

The features of the T-Z trace (Figure 34a) are easily 
described by the G-Z model. The product of signed strength and magnitude 
gives the total temperature change from the surface temperature to a 
given depth for a step-structure temperature profile. The temperature 
at any depth could then be found if the sea surface temperature was 
known* Thus, with magnitude, signed strength and sea surface temperature, 
dny layered or stepped thermal profile could be reconstructed. 
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Figure 32. Two-Layer Model. 
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Figure 34. Step Structure Model. 
65 



2 . Continuous Curving Gradients 



Since finite values are used for both depth and temperature, a 
truly continuous curve cannot be drawn. The spacing by the CalComp 
plotter of every hundredth of an inch can cause variations of this 
size. Computed values of depth, temperature and temperature gradient 
were provided as inputs to a CalComp plotter at every hundredth of an 
inch to create a quasi-continuous curve, at least to a resolution of the 
plotter. It is this quasi-continuous field which accounts for the 0.01 
inch steps encountered in curves drawn by the plotter for any truly con- 
tinuous curve. This effect is also present whenever the plotter is 
drawing any line not directly in the X or Y axis directions. Accordingly, 
any steps shown of this size could be considered as possible adjustments 
made by the plotter or data to adjust to the hundredth of an inch grid. 

Two examples of continuous gradients are given to show the 
resulting distributions of gradient with respect to depth, 
a. Gaussian 

An example of a Gaussian profile (Figure 35) shows a normal 
distribution of temperature with depth for the function, like 

Z 1 2 

20 r ” 2 ^ 

= 80 - ^ J e dx (8) 

-00 

\ 

where: x = a dummy variable 

80 = sea surface temperature. 

The G-Z curve highlights the changes in the T-Z curve. The increase in 
gradient strength is very clear, reaching a maximum 200 feet, whereas 
the inflection point of the T-Z plot is difficult to find by visual 
inspection. The summation of signed strength values at every foot of 
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Figure 35. Gaussian Model. 
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depth from the surface down gives a total temperature change which is 
analogous to the product of magnitude and signed strength with all 
magnitudes set equal to one. With this summation and a given sea surface 
temperature (SST) , a continuous temperature profile can be reconstructed, 
b. Exponential 

Using an exponential expression for the temperature distri- 



bution 



T(Z) = 60 



20,-ttz/iooo 



(9) 



where : 



60 = SST, 



the T-Z (Figure 36a) and G-Z curves (Figure 

of each other. The derivative of T(Z) with 

dT(Z) _ 
dZ 50 



36b) appear as compliments 
respect to Z gives : 

( 10 ) 



which is the equation of the G-Z plot. 

3 . Superpositions 

In the example of layers superimposed on an exponential model 
(Figure 37), the conditions for the exponential (9) hold for the depth 
intervals: surface to 40 ft, 220 to 340 ft and 460 to 800 ft, with 

various layers at the remaining depths. Isothermal layers exist below 
the weak surface exponential gradient layer and also in the 140 to 180 ft 
depth layer. From the T-Z plot, the layers at 340 and 400 ft are not as 
evident as they are in the G-Z plot. 

4 . Invers ions 

A model with two inversions and several linear steps (Figure 38) 
shows the positive gradient regions clearly. The second deeper inversion 
is clearly a weak gradient, but it is still discernible in Figure 38b. 
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Figure 36. Exponential Model. 
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Figure 37. Exponential Model with Layers Superimposed. 
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Figure 38. Multi-Step Model With Two Inversions 
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With inversions, summing the signed strength values for each foot gives 
the temperature change from the surface and adding this value to SST 
yields the actual temperature at any depth. 

B. OBSERVATIONAL CASES 

Several different oceanic regions were selected to establish fea- 
tures of this model. The oceanic regions explored were along the Atlantic 
Coast^in the Andaman Sea, Eastern North Pacific, and the Gulf Stream. 

1 . Atlantic Coastal Region 

Stations in Massachusetts Bay occupied during December 1969 
were used. The stations were within a radius of 13 miles. Water depths 
were between 280 and 420 feet. Lack of strong gradients in all seven 
traces shows that the water was probably well mixed (Figures 39 through 
42). A positively signed strength in the upper 300 feet was common in 
the traces. Inversions (Figure 40) are quite clear in both the T-Z and 
G-Z plots. A portion of the actual XBT trace is shown for three of the 
stations. In these three cases, there is excellent visible correlation 
between the actual XBT trace and the reconstituted T-Z plots which shows 
that MICROXBT can resolve very small temperature features. In no case 
was there an isothermal layer greater than 140 feet thick. 

2 . Andaman Sea 

In the Andaman Sea, near the Straits of Maalacca, a series of 
15 XBT stations from the 18 to 20 July 1969 period were used. The 
stations were within a radius of 12 miles. Water depths varied between 
290 and 380 feet. This oceanic region could be considered somewhat 
similar to the Atlantic Coastal region above. 
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jure 39. T-Z and G-Z Plots for Atlantic Coastal Region 6 and 7 December 1969. 
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Figure 40. T-Z and G-Z Plots for Atlantic Coastal Region 7 December 1969. 
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;ure 41. T-Z and G-Z plots for Atlantic Coastal Region 9 and 10 December 1969. 
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Figure 42. 



T-Z and G-Z Plots for Atlantic Coastal Region 



11 December 1969. 
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The T-Z plots (Figures 43 through 50) contain a mixed layer to 
200 or more feet except for one observation at 1200Z on the 19th (Figure 
48b). The mixed layer transients are marked with spikes in the G-Z plots. 
Figure 50 shows transients with the greatest magnitude in the upper 100 
feet. The seasonal thermocline is characterized by two steps (Figure 
44, for example). The G-Z plots in most cases show this second step, 
averaging at a depth about 300 feet, with good resolution. Some strong 
gradients are found in the 200 to 300 feet depth region in every G-Z 
plot except in Figure 48b. This figure appears to indicate flow over a 
hump or have erroneous depth values. The fall-rate of the XBT probe, 
appearing erroneously in the T-Z plot, indicates a possible malfunction — 
of the probe or the recorder for this one cast. 

3 . Pacific Oceanic Region 

A series of XBT traces, taken over a week period in October 1968, 
in the deep waters (over 100 nm off the West Coast near San Diego, 
California) were used for this observation. A composite T-Z plot of 
the actual XBT temperature-depth traces is shown in Figure 51. This 
region is sub-tropical oceanic [Tully, 1964]. Plots of T-Z and G-Z 
(Figures 52 through 65) are presented from the surface to 900 feet. The^ — 
seasonal thermocline is evident in all plots. The strong thermocline is 
clearly shown. Some transients are notable. Few inversion layers were 
present, and, generally, they were found only in the near-surface layer. 
The G-Z plot is characterized by strong gradients in the seasonal thermo- 
cline and by weak gradients in the main thermocline, each separate layer 
separated by isothermal waters from one to forty feet thick. 

The mixed layer contained numerous temperature fluctuations in 
most cases. Some inversions, in particular the one at 0600Z on 14 October 
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Figure 43 . T-Z and G-Z Plots for Andaman Sea Region for 18 July 1969. 
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ure 44. T-Z and G-Z Plots for Andaman Sea Region for 18 July 1969. 
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Figure 45. T-Z and G-Z Plots for Andaman Sea Region for 18 and 19 July 1969. 
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lire 46. T-Z and G-Z Plots for Andaman Sea Region for 19 July 1969. 
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Figure 50. T-Z and G-Z plots for Andaman Sea Region for 20 July 1969. 
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Figure 51. Composite T-Z Plot from Actual XBT Traces Down to 900 feet 
for Region 100 nm off San Diego. Data are given In time 
sequence from 10 through 16 October 1968. 
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Figure 56. T-Z and G-Z Plots for Eastern North Pacific for 1800Z 12 October 1968. 
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Figure 58. T-Z and G-Z Plots for Eastern North Pacific for 0600Z 13 October 1968. 
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igure 60. T-Z and G-Z Plots for Eastern North Pacific for 1800Z 13 October 1968. 
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Lgure 62. T-Z and G-Z Plots for Eastern North Pacific for 1200Z 14 October 1968. 
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1968 at 120 feet was clear on both the bathythermogram and the T-Z plot. 
The steps are quite irregular on the T-Z plot but not as irregular as on 
the G-Z plot. 

4 . Gulf Stream Region 

Thirteen stations, during the period 11-13 August 1969, were 
analyzed for the Gulf Stream region (Figures 66 through 80). The ship*s 
track crossed the Gulf Stream boundary several times. 

a. Interior of Gulf Stream 

For the interior of the Gulf Stream (Figures 70 and 71), 
the mixed layer is distinct and a typical seasonal and main thermocline 
are shown. A "boot” in the 100 foot layer is clearly shown by positive 
gradient values in several of the figures and the actual XBT traces. 

b. Frontal Region 

Strong gradients are the rule for this region and they are 
quite distinct in both the T-Z and G-Z plots. The contrast between 
the frontal and interior regions of the Gulf Stream are quite apparent. 
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Figure 66. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T-Z Trace Inset. 
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Figure 68. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream, Actual XBT T-Z 
Trace Inset. 
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Figure 69. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T-Z 
Trace Inset, 
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Figure 70. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T-Z 
Trace Inset. 
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Figure 71. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T-Z 
Trace Inset . 
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Figure 72. T-Z artd^Z Profile for Transects across Northern Boundary of the Gulf Stream. Actual XBT T~Z Trace Inset. 
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e 73. T-Z and G-Z Profile for Transects' across Northern Boundary of the Gulf Stream. Actual XBT T-Z Trace Inset 
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Figure 76. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 77. T>Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 78-:- T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 79. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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Figure 80. T-Z and G-Z Profile for Transects across Northern Boundary of the Gulf Stream. 
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VI. CONCLUSIONS 



The results give a measure of the microthermal structure obtained 
with an XBT by the NDM. 

A. PHYSICAL SIGNIFICANCE 

Microthermal structure statistics have not been adequately explored. 
This thesis has proposed an easy method to observe the fine-scale thermal 
structure. Further, this proposed objective analysis scheme does not use 
the T-Z values but rather the G-Z values. It is the G-Z distribution 
which is important. The T-Z can always be reconstructed if one additional 
parameter, the SST, is known. 

1 . Internal Waves - Vaisala Frequency 

The Brunt-Vaisala frequency is the maximum frequency of oscil- 
lation of free internal waves. It may be calculated from the density 
gradient, or from the parameters which define density gradient. In 
isohaline water, the NDM gives a direct measure of density (or stability). 
In such waters it could be used to set an upper limit on the maximum 
frequency of oscillation of internal waves. Further, a time series 
study would clearly illustrate the amplitudes likely to be encountered 
since the sheets indicated on the G-Z diagram are readily identified. 

A correlation study of G-Z, Brunt -Vaisala versus depth diagrams would 
be useful. 

2 . Oceanic Fronts 

This difficult region to analyze microscopically could possible 
be better explored by using the NDM for the microthermal processes. The 
NDM explores the continuously changing conditions in frontal regions. 
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Frontal regions are usually indicated by extremely complex G-Z plots 
which illustrate the extreme gradients, both positive and negative, 
encountered. Based on the admittedly small number of cases examined, 
there appears to be a tendency of the layers to maintain their character- 
istics and slide by, over or under each other rather than mix. 

3 . Sound Propagation 

The effect of microthermal structure on sound propagation 
require better statistical knowledge of the oceans. The NDM gives 
emphasis to the strongest gradient regions. 

The microthermal effects on index or refraction variations 
(density fluctuations) must be determined. Using the NDM with salinity 
variations (from an STD cast) and the NDM with sound velocity variations 
(from an SVTP cast) taken from the two instruments simultaneously, cross- 
correlation techniques could shed some light on the unresolved micro- 
processes in the ocean. In particular, what are the exact physical and 
theoretical correlations among the salinity, temperature, and sound 
velocity. 



B. SUGGESTED STUDIES 



Possible studies for the 
developing this thesis. 

1 . Mathematical Models 
Various mathematical 
because each method had some 
their utility. 

a. Analysis Methods 
The microthermal 
as a ’’roughness” factor, the 



future have suggested themselves while 



models attempted could be explored further 
merit. Suggested models are limited by 



structure recorded by XBTs can be looked at 
noise level of the water structure. Further 
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work is required to determine such a factor for practical operational use. 
This roughness factor should be computed as the XBT trace is being 
analyzed at FNWC by a program like DIGITBT. This factor could be one 
number to describe objectively the microthermal structure variations on 
the G-Z trace. 

Variations in the NDM should be used with the variables that 
have been arbitrarily set. For example, the depth interval of AT/ft 
at one foot should be varied to find the effects a larger or smaller 
depth interval would have on the G-Z profiles. As the steps become 
smaller, the exact derivative is approached. Larger steps should smooth 
the spikes shown in the G-Z figures. In some cases, this smoothed G-Z 
profile might be described with curve fitting techniques as suggested by 
Boston [1966] and Grosfils [1968]. 

MICROXBT^s digitizing interval of 0.01 inches should be 
changed to 0.005 inches and the program modified to use this interval. 

This would yield more detailed thermal structure than MICROXBT. Obviously, 
one cannot make the interval finer than the noise level of the sensor or 
recording instrument . 

b. Physical Processes 

The G-Z profiles from the NDM emphasize the changing physical 
features which are obscure in the T-Z profiles. By using G-Z profiles 
for parameters other than temperature, clarity in determining distri- 
butions of chemical parameters would be similarly emphasized. Limitations 
in use of the NDM on other chemical properties should be further 
investigated . 

Water mass identification is possible with G-Z profiles. 
Additional work is required in this area. 
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2 . Experimental Work 



A microthermal sheet probe could be designed for investigating the 
chemical and physical processes inside the sheet and inside the layers. 

A comparison of actual physical property differences should be made. 

A modified expendable temperature gradient (XTG) meter could be 
designed using two thermistors to obtain the temperature differences in 
the probe. The standard XBT system equipment already installed could be 
changed to give analog G-Z profiles to the operational user without any 
de lay . 

a. Time Series 

The XTG meter could be used to construct G-Z variability at 
one location by conducting a time series analysis. This needs to be 
done in order to obtain a roughness factor, an objective microstructure 
measuring parameter . 

b. Transect 

An XTG meter could be used in transects across various water 
masses to obtain the definite characteristics for the water masses. The 
XTG could be used to follow interesting features on the G-Z plots. More 
important, however, is that the special characteristics of microthermal 
structure could easily be determined. 

3 . Handling of Data 

When the microthermal structure roughness factor is determined, 
this factor should be made an integral part of the records for that XBT 
trace. This information should be included as part of the ID data. 

Aboard Navy ships the AN/UNQ-7 tape recorder could be used to 
obtain magnetic tape records of the analog XBT or XTG signals. These 
tapes could be submitted with the traces to FNWC or NFS for further analys 
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4 . Naval Applications 



Results from microthermal studies show that high frequency sonar 
performance as influenced by the microstructure receive the greatest 
microstructure effects. More knowledge of the microstructure is needed 
before the torpedo and mine-hunting sonar performance could be improved. 
With a measure of the microstructure, the performance of a sonar system 
could be predetermined. 

Tens of thousands of XBT traces are on file at FNWC . New data 
are arriving daily. This thesis is an attempt to show how this data can 
be used both operationally and scientifically. Obviously it represents 
only a beginning in the way of microstructure analysis and operational 
utilization. 
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APPENDIX A 



MICROXBT FORTRAN IV PROGRAM FOR CDC 6500 
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♦** COMPARE TRACE RECORDS AND MERGE IF COMPATIBLE. TRACES ARE MERGED A195 
♦♦♦ IF FOUND TO BE WITHIN THREE INCREMENTS OF EACH OTHER. A196 
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1021 F0RMAT(X,R4,X,6R1,X,4R1,X,5R1,X,6R1) ft055 

IIS(IS)=IT(1) $ IIS(IS+1)=IT(2) $ IIS(IS+2)=IT(3) A519 

IS=IS+3 $ ICT=ICT+1 A520 
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PRINT 1031, ICT A554 

1031 FORMAT(23H SATISFACTORY TRACES ,14) A065 

STOP A555 

END A556 
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13 IP(H)=IX *100 S IPCM+ 1 )=IY *100 A 633 

14 IF(IP(M) .EQ.IP(M- 2 ) .ANO.IP(M+l).EQ.IP(M-in GO TO 6 A 634 

IF( IP(M) .GE. IH( 40 ) ) GO TO 7 A 635 

M=H +2 A 636 

GO TO 5 A 637 
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APPENDIX C 



JOB CONTROL LANGUAGE (JCL) FOR READING FNWC TAPES 
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//(Job name) J0B (Acct. no., Job no., Section), '(Your name)' Note 
// EXEC F0RTCLGP,REGION=1OOK,TIME.G0=2 Note 
//F0RT.SYSIN DD * 

F0RTRAN Source Deck 

/* Note 
//G0.FTO6FOO1 DD SYS0UT=A,SPACE (CYL, 6) Note 
//G0.FTO4FOO1 DD UNIT=24OO-1,V0L=SER=NPSXXX,IABEL=(,NL) , Note 
// DCB=(DEN=1,RECFM=F,LRECL=136,BLKSIZE=136,TRTCH=ET) , " 
// DISP=(0LD,KEEP) Note 
//G0.FTO4FOO2 DD AFT=FT04F001 ,V0L=SER=NPSXXX , IABEL(2 ,NL) , Note 
// DCB=(DEN=1,RECFM=F,LRECL=136,BLKSIZE=136,TRTCH=ET), " 
// DISP=(0LD,KEEP) Note 
//G0.SYSIN DD * ' , Note 



DATA Cards, if any 

/* (Orange Card) 

Notes 

1. Standard J0B, green card, for program submission. 

2. If not plotting or card punching, the F0RTCLGP can be changed 
to F0RTCLG. 

3. Must be a white card. 

4. Optional card for additional storage for output. 

5. The XXX following NFS is the NFS tape number. These three cards 
are for reading the first file on the input tape. These cards 
assign input unit 4 for the source programs READ statements. 
Further note that the record length is 136 characters, all of 
which may be read. This is the standard record length from the 
CDC 3600. 

6. For a third file to be read, these three cards must be repeated 
with F002 changed to F003 and LABEL(2,NL) changed to IABEL(3,NL) 
and so forth for additional files. Also the UNIT=2400-1 should 
be changed to AFT=FT04F001 for all files after the first file. 

If only one file is to be read, then these three cards are not 
required . 

7. Not required if there are no input data cards. 



1 

2 

3 

4 

5 

5 

5 

6 

6 

6 

7 
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APPENDIX D 



JOB CONTROL LANGUAGE (JCL) FOR DUMPING FNWC TAPE DATA 



14A 



//(Job name) 


J0B (Acct. no., Job no., Section), ’(Your name)' 


Note 


1 


//DUMPTAPE 


EXEC PGM=IEBPTPCH 






//SYS PRINT 


DD SYS0UT=A 






//SYSUTl 


DD UNIT=2400- 1 , DISB=0LD , LABE L= ( , NL) , V0LUME=SER=NPSXX , 


Note 


2 


// 


DCB= (DE N= 1 , RECFM=F , LREC L=136,BLKSIZE=136, TRTCH=ET ) 


Note 


2 


//SYSUT2 


DD SYS0UT=A 






//SYS IN 


DD *,DCB=BLKSIZE=80 








PRINT MAXFLDS=1,ST0PAFT=1OOO 


Note 


3 




RECORD FIELDS=(120) 


Note 


4 



/* (Orange Card) 



Notes 



1. Standard J0B, green card, for program submission. 

2. If more than one file is on the tape then separate two cards for 
each file identified by file number by the file number in the LABEL 
(X,NL) where S is the file number, which is implied 1 when omitted. 

3. The ST0PAFT limits the number of lines of printed output, in this 
example 1000 lines is used before printing is terminated. 

4. The field is recorded at the maximum amount that this program will 
dump which is 120 characters. If there are any characters between 
120 and 133 in the records then this routine will not dump those 
characters . 
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